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Recent inelastic neutron scattering experiments in Celna and CePd2Si2 single crystals, measured 
spin wave excitations at low temperatures. These two heavy fermion compounds exhibit anti- 
ferromagnetic long-range order, but a strong competition between the Ruderman-Kittel-Kasuya- 
Yosida(RKKY) interaction and Kondo effect is evidenced by their nearly equal Neel and Kondo 
temperatures. Our aim is to show how magnons such as measured in the antiferromagnetic phase of 
these Ce compounds, can be described with a microscopic Heisenberg-Kondo model as introduced 
by J.R.Iglesias, C.Lacroix and B.Coqblin, used before for studies of the non-magnetic phase. The 
model includes the correlated Ce-4/ electrons hybridized with the conduction band, and we consider 
competing RKKY (Heisenberg-like Jh) and Kondo (Jk) antiferromagnetic couplings. Carrying on 
a series of unitary transformations, we perturbatively derive a second-order effective Hamiltonian 
which, projected onto the antiferromagnetic electron ground state, describes the spin wave excita- 
tions, renormalized by their interaction with correlated itinerant electrons. We numerically study 
how the different parameters of the model influence the renormalization of the magnons, yielding 
useful information for the analysis of inelastic neutron scattering experiments in antiferromagnetic 
heavy fermion compounds. We also compare our results with available experimental data, finding 
good agreement with the spin wave measurements in cubic Celna. 



PACS numbers: 71.27.-|-a,75.30.Ds 



I. INTRODUCTION 



The description of heavy fermion compounds is chal- 
lenging due to the rich variety of phase diagrams they 
present, and the anomalous physical properties which 
may be found. Among them, appear the Ce and U 
compounds which exhibit long-range antiferromagnetism 
(AF) at low temperatures (for example, among antiferro- 
magnetic Ce compounds, CeRh2Si2 exhibits the highest 
ordering temperature T/v = 36K, with local magnetic 
ordered moments of 1.34 — 1.42/ib per Ce, i.e. rela- 
tively large, compared to the full Ce^+ free-ion value: 
2.54 ^ib)'^ Depending on the particular compoundfSi the 
antiferromagnetism takes different forms (magnitude of 
the local moments varies widely: e.g. 0.001 /is as 
in CcRu2Si2 or 0.02 ^xb in UPts, to 1.55 ^xb as in 
UCu5; as does the spin configuration: with three-, two- 
or one-dimensional AF structures observed). Antifer- 
romagnetism may also appear competing or coexisting 
with superconductivity^^!^ for which spin fluctuation- 
mediated pairing mechanisms are explored;^^ Non-Fermi 
liquid behaviour may appear;^ and quantum criticality 
has become a subject of intensive study in these com- 
pounds, both experimentally and theoreticallyj^i^i^iiSiiii 
The crossover from the antiferromagnetic state to the 
non-magnetic heavy fermion state, which can be tuned by 
pressure, doping or magnetic field, is one of the most in- 
teresting problems in strongly correlated /—compounds. 

The physical properties of these compounds are deter- 
mined by the strongly correlated /—electrons present (4/ 
in Ce; 5/ in U) and their hybridization with the conduc- 
tion band. The RKKY indirect exchange interaction be- 



tween /— local magnetic moments, favouring the estab- 
lishment of long-range magnetic order, competes with the 
screening of these moments by the conduction electrons, 
described by the Kondo effect This competition is the 
subject of the Doniach diagram, ^'^ which compares the 
variation of the Neel and Kondo-impurity temperatures 
with increasing antiferromagnetic intrasite exchange cou- 
pling Jk , between local /-moments and conduction elec- 
tron spins. Compounds with similar magnetic order- 
ing temperature T/v and Kondo temperature Tr-, the 
temperature below which magnetic susceptibility satu- 
rates indicating coherent Kondo-singlet formation, are 
ideally suited to the study of this RKKY-Kondo compe- 
tition. In this regard, experiments on CeM2Sn2 (M=Ni, 
Ir, Cu, Rh, Pd, and Pt: Ad or hd transition metals)i^ 
and CeX2Si2 (X=Au, Pd, Rh, Ru)— were undertaken, 
indicating that departures between theory and experi- 
ments resulted from the use of Kondo impurity relations. 
The Kondo-lattice model, instead, consisting of a lattice 
of local magnetic moments coexisting with a conduction 
band, has proved appropriate for the description of many 
4 — / and 5 — / materials, in particular most Ce (or Yb) 
compounds, respectively corresponding to a configura- 
tion close to 4/^ (or 4/"'^^), where one 4/ electron (or 
hole) interacts with the conduction elect ronsJ^ii^iii!^ 
In 1997 a revisited Doniach diagram was introduced, in- 
cluding short-range antiferromagnetic correlations in the 
Kondo lattice, in order to improve the description and, 
in particular, to account for the observed pressure depen- 
dence of T/f in CeRh2Si2»^ The situation is more com- 
plex in Uranium compounds, where U has a 5/" configu- 
ration with n=2 or 3, since the 5/ electrons are much less 
localized than the 4/ electrons of rare earths. Regarding 
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spin dynamics, it is not clear that Ce and U compounds 
are intrinsically similar*^ In the following we will focus 
on Ce-compounds, except otherwise specifically stated. 

Experimentally, while the magnetic response due to 
Kondo spin fluctuations in the paramagnetic state of 
heavy fermions is well studied, relatively little is known 
about the nature of the magnetic excitations in the or- 
dered phase of Kondo lattices,^ on which our present 
study will focus. Being still an unsolved problem how 
to describe on equal terms both the Kondo effect and 
antiferromagnetism,JS here we will focus on systems with 
relatively large local moments and study them deep in- 
side the antiferromagnetic phase: far from the antiferro- 
magnetic quantum critical point, where spin fluctuations 
would become more relevant. 

A few years ago CePd2Si2^' single crystals were stud- 
ied with inelastic neutron scattering: below the antifer- 
romagnetic ordering temperature strongly dispersive spin 
wave excitations were found, with an anisotropic damp- 
ing, which coexisted with the Kondo-type spin fluctua- 
tions also present above Tjy. At T = 1.5K these spin 
waves were measured along various BZ paths: they were 
found to present an energy gap of 0.83 meV and to ex- 
tend up to almost 3.5 meV. CePd2Si2 has a bcc tetrag- 
onal structure, and its antiferromagnetic ground state is 
characterized by propagation vector: q — (1/2,1/2,0), 
with ordered moments: m = 0.66/iB, T^v = 8.5K and 
Tk = lOK, and linear electronic specific heat coef- 
ficient 7 — 250mJ/molK^. Under pressure applica- 
tion, at 28.6 kbar the system undergoes a transition 
into a superconducting phase with critical temperature 
of 430 mKi2ii2^ More recently, inelastic neutron stud- 
ies of Celna single crystals were performed, with simi- 
lar resultsi^ Well defined spin wave excitations with a 
bandwidth of 2meV and a gap of 1.28 meV were found 
in the antiferromagnetic phase^^ coexisting with Kondo- 
type spin fluctuations and crystal-field excitations which 
also appeared above T/v = lOK = Tk- Celns crystal- 
lizes in a cubic (fee) structure, with an antiferromagnetic 
structure characterized by magnetic propagation vector 
q = (1/2,1/2,1/2), with ordered moments: m = 0.5/ib 
and 7 = 130mJ/molK^. Under application of pressure, 
at 26.5 kbar the system undergoes a transition into a 
superconducting phase with critical temperature of 200 
mK.^ 

In next section, we will briefly introduce the micro- 
scopic Heisenberg-Kondo model proposed by J.R.Iglesias, 
C.Lacroix and B.Coqblin,— to study the non-magnetic 
phase of heavy fermion AF compounds, to which we 
shall add conduction electron correlations. We will then 
present our calculation for the renormalization of spin 
wave excitations due to their interaction with the corre- 
lated conduction electrons (the Appendix complements 
this section). In Section IIIIl we will discuss the re- 
sults of our study, show how the different parameters of 
the model influence the renormalization of the magnons, 
and compare our results with the available experimental 
dataf22i^ In Section lTVl we summarize and point out that 



the present work should yield useful information for the 
analysis and prediction of inelastic neutron scattering ex- 
periments in heavy fermion AF compounds, as CeRh2Si2. 



II. MICROSCOPIC MODEL, AND 
PERTURBATIVE APPROACH. 

In order to describe the Ce-heavy fermion systems ex- 
hibiting antiferromagnetic long-range order we have used 
the microscopic model which has been proposed by Co- 
qblin et ali^ to describe the competition between the 
Kondo effect and the RKKY interaction, in compounds 
where departures from the original Doniach picture^ ap- 
pear. In principle, both the RKKY magnetic coupling 
and the Kondo effect can be obtained from the Kondo 
intrasite-exchange term, but when dealing with approxi- 
mations it is difficult to insure that both effects are taken 
into account if an explicit intersite exchange (as the effec- 
tive RKKY interaction or, depending on the system, also 
the direct exchange) is not included in the Hamiltonian.^" 
The model-^consists of a Kondo lattice, featuring local 
magnetic moments coupled both to conduction electrons, 
by a Kondo-type interaction , and among themselves 
by an antiferromagnetic RKKY-type exchange > 0. 
The moments are assumed to order below the Neel tem- 
perature T^r. The Hubbard-correlated conduction elec- 
trons occupy a non-degenerate band. Therefore, the 
model may be represented in standard notation by the 
following Heisenberg-Kondo Hamiltonian: 

H ~ ^bajid + ^Kondo + Hneis (1) 
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H Kondo — Jk \ ■ + ' Sj" 

\i<^A jet3 



(3) 
(4) 



Here I (j) indicates that the lattice site index j runs 
over the z nearest neighbours of site I. The itinerant elec- 
tron spin is Sj , while A, B label the two interpenetrating 
sublattices of local moments, >5';g_4 and S^^q , with op- 
posite moment direction. The Kondo exchange Jk could 
in principle have either sign (though for heavy fermion 
compounds, it would be antiferromagnetic: Jk > 0). In 
the following we shall distinguish in the Kondo term the 
longitudinal H^^ from the transverse Hj^ contributions. 
By taking as positive z direction the direction of >S';g_^ , 
they are defined as: 



(5) 
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leA jet3 
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We diagonalize the Heisenberg term, by represent- 
ing the local moments operators in the Holstein- Pri- 
makoff approximation, which is appropriate at tempera- 
tures much lower than Tjv . Namely, we take 



-5*; - {S^,Si ,Si) 

^ {V2Sbi,V2SblS-blbi) 
Sj - {V2Sb],V2Sb^,-{S-b]b, 



(6) 



where 6^'' and b'^J^ are the bosonic operators in sublat- 
tices A and B, respectively. Going to reciprocal space in 
the reduced Brillouin zone (RBZ) and introducing the AF 
spin wave operators {agjOg} by the Bogolyubov trans- 
formation: 

at = fotch(^,)-6_,Sh(^,) 

b^Ch{^,)-bl^Sh{^,) (7) 

Th{2i}g)^--J2cos{qAij) (8) 



where 



diagonalizes Hueis- In the absence of the interaction 
with the fermions, i.e. in the limit of vanishing Jk, the 
frequency of the bare AF spin waves would be {z is the 
number of nearest neighbors and Aij is the vector joining 
two n.n. sites) 



HHeis = ( "9"? + I 



Ea,, cos {qAij] 



(9) 



We further assume that the band electrons have devel- 
oped an AF spin order, due to the cooperating effects of 
the Hubbard correlation U and of the interactions with 
the AF-ordered local moments, which provide a staggered 
field. 

To diagonalize Hband in the AF state, we will use a re- 
formulation of Gutzwiller's variational approach for the 
description of antiferromagnetism in narrow bands due 
to Spalek et al.— This approach allows to connect the 
standard (mean-field) Slater band-insulator to the local- 
ized Mott antiferromagnetic insulator. One expresses the 
correlation-induced bandwidth reduction in the param- 
agnetic (PM) state by a Gutzwiller-type factor $(71,77) 
depending on the band filling 77 and on the probability of 

double occupancy 77 — N^^ ^77;|77;|^ . The correlated 

band energies are written as = $(71, 7/)e^, where 
are the uncorrelated band energies and: 



$(77, 77) = 1 - 



1 i 



(10) 



The U- depending optimal value of 77 is found at zero 
temperature by minimizing the PM energy 



E 



PM 



^$(77,77)4 (77 



kal 



NUt] 



at given 7i and U . 

Assuming that the electrons have an AF ground state 
of Neel-type one adopts the standard Slater formalism, 
only with the PM energies renormalized according to 
EqdOl 

We introduce the fermion operators for this AF Slater- 
type state IttfeCT'/Sfclj^l by the transformation 



-^k+Q,a 



PL cosCfc^ -aL sinCfca 
/3L sin Cfc<T + ala cos Cfc^ 



(11) 



where Q =(1/2, 1/2, 1/2) in units of 271 /a is the wavevec- 
tor characterizing the AF magnetic state and k is 
a wavevector belonging to the reduced Brillouin zone 
(RBZ) defined by |k| < |Q/2|. The diagonalization con- 
dition for a lattice with an inversion center yields: 



tan (2Cfc<T) 



fT tan {2C,k 



(12) 



where the AF order parameter (sf 



Is) of the itiner- 



ant electrons (staggered magnetization) , or band electron 
polarization, is given by: 

1 -^r^ 



2N ^ 

keRBZ,a 



'■ka/ 



^ka 



sin 2Ck 



Therefore, the diagonal bare electron Hamiltonian reads: 



ka k<T 



keRBZ,a 



(13) 

where, assuming that the lattice has a center of in- 
version, the bare electron eigenenergies (actually spin- 
independent) read {x = a,f3): 



1 



<^^el + {U{s))\ (14) 



where 5xa denotes the Kronecker delta. Notice that 
a is the lower subband. The subband filling factors 

7f _ \ = [exp [El, - /ksT + l]-\ho de- 



pend on (s) through E'^ 



ka ' 



A. Renormalization of the bare subsystems by part 

of Hi, 



By expressing the longitudinal Kondo term in terms of 
the AF bare electron operators | and the bare 

bosonic operators |agt^|, it can be rewritten and decom- 
posed into contributions which are respectively diagonal 
in the bare bosons, i.e. H^'^\ or non-diagonal, i.e. E. 



Explicitly, the diagonal term reads: 



^Sj2[^in{2Ck){nt,--n>, 
-a cos (24.) [alp,, + f3la,,)] [1 - V,] 



(15) 
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where the bosonic operator Vq reads: 



NS 



iVS* 



(16) 



The term non-diagonal in the bare bosons, wiU be 
discussed in next section. 

Now, we treat H^'^'^ in mean-field approximation 
(MFA), obtaining three kinds of terms. The first one, 
obtained averaging over the bare electrons reads: 



H 



z{dl} 
K 



~Jk (s) (ala^ + al^a) Ch (2^,) 



(17) 



Jk (s) ^ [ala'ig + a^a^A Sh (2i9q) + const 



and renormalizes the spin wave frequencies, therefore 
modifying the diagonalization condition (Eq. ^) for the 
AF Heisenberg Hamiltonian as: 



Th {2dq 



(18) 



Notice that, in order to minimize the mean-energy of the 
Kondo term, —J^ (*) is a non- negative quantity, irre- 
spective of the sign of Jj^ ■ 

The MFA-redefined Heisenberg term now reads: 



HHe^s=Y.^^<ly'W+2) ^^^^ 
q ^ ^ 

with the MFA-renormalized frequency given by: 
hn^ EE {zJhS+ I Jk {s) I) . 1 - 



Ea,, cos(gAy) 



Y [z+\Jk{s)\/JhS 

(20) 

Averaging over the bare bosons yields two terms. To 
write them it is convenient to define the local moment 
reduction factor as follows: 



-^[(aX)Ch(2^,)-HSh2(^,)] (21) 



where the mean number of MFA bosons is (^a^gaqj = 

[exp (ftrig/fesT) — 1] ^. Then the (in general 

temperature-dependent) expectation value of the local 
moment is {S^) — S — Afsw 

The second contribution, which reads 



(22) 



renormalizes the electron eigenenergies. The diagonal- 
ization condition for the AF electron Hamiltonian now 



tan(2Cfc<T) = -CT 



U {s) - Jk jS^) /2 



(23) 



Consequently, the AF order parameter (s) of the itinerant 
electrons is given by: 



2N 



^ [U is) - Jk (S^-) /2] - 
k. ^^^el + [U{s)-JK{S^)/2f 



and the MFA eigenenergies by {x — a,P): 
E 

z(d3) 



(24) 



L = \un + (1 - 26.^) ^<^^el + [U{s)~JK{S^) I2f 

(25) 

The third contribution H^j^""^' provides a hybridization 
between the a, (3 electron states 



The transformation 



"L = ('OS £,ka+Bl^ sin ^ka 

PL = Bl^COS^ka- sin ^kc 



where: 

tan (2^,7) = cr 



-Jx(5^)cos (2Cfc) 



E^-E^- Jk{S-) sin (2^) 



(26) 



(27) 



(28) 



allows to eliminate H^'^^^ yielding the MFA hybridized 
electron Hamiltonian in diagonal form as: 



ttAF 
band 



(29) 



The hybridized band energies are {X — A, B) 



cX _ 

Oi. — 



1 



E^ + E'^)-i26xA-l 



(30) 



where 

Tk = 



Et E^k 



1 - 



Jk{S') sin(2Cfc) 
E2-ES 



Jk {S') cos (2Cfe) 



E^ - E^, 



(31) 



k,cr 



is the difference of energy at wavector k between the AF 
subbands. Notice that the A band is the lower one. Tak- 
ing the paramagnetic band filling per site n < 1 only this 
band will be non-empty. 
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B. Perturbative derivation of the effective magnon 
Hamiltonian. 



At this stage, we will describe the perturbative treat- 
ment performed to derive the effective second-order 
Hamiltonian for magnons. We begin by rearranging the 
Hamiltonian terms, splitting it into a basic part Ha, in- 
cluding the MFA redefined "bare" magnon and electron 
Hamiltonians of previous section, plus a "perturbation" 
J, which includes the full transverse Kondo coupling term 
(/-'- = H^) and the scattering (i.e. non-diagonal in bo- 
son operators) part of the longitudinal Kondo coupling 
(/^) not considered so far. 

Explicitly, we have H = Hq + I = Hq + + I-^ where 



Ho = HHe 

reads: 



H 



AF 
band' 



The longitudinal perturbation 



pqr,a 



(^rcr^p+l — q.cr ^rcr^p+r-q,a) ^Pq 

Jk 



N 



Cr (1 — Spq) L^^p-^r-q 



X (A|!^Bp_|_^_q + Bl^Ap^j,_ij^^) Gpq (32) 
The spin wave operators are contained in 

Gpq = ala^Ch{^q)Ch{^p) + alqa^pSh{^q)Sh{^p) 
+alalgSh i^q) Ch (dp) + a_pa^Ch (dq) Sh (i9p) 

while coefficients M^^j^^_^ and L^p^^-q defined in 
the Appendix. 

The transverse perturbation is written as: 

^ J^yi^ { KCh i^q) + a_^Sh i^q)] (s+ + 

+ [aJSh i^q) + a_^Ch i^q)] [stqs + SZq,) } (33) 

The spin operators are expressed through the electron 
operators as: 



^qA ^ "qB 



k+q.-a^BB i^' ?) 



l.flyBl B 

-l\[^T.''BlA+q,-.C'nAik.q) 



(34) 

where A,t = ± and the coefficients C^'^(fc,g) with 
X,Y = A, B are defined in the Appendix. 



In the previous section, we have included at MFA level 
the effects of the diagonal part of the longitudinal Kondo 
coupling on the bosonic states (Eqs. [T9ll20l ) and on the 
AF electronic states (Eqs. [5011511 ). We will now deter- 
mine the magnon renormalization effects derived from 
the transverse and non-diagonal parts of the longitudinal 
Kondo coupling, which we are considering as the pertur- 
bation / on the MFA state. 

The effect of the perturbation will be taken into 
account through a Frohlich-type of truncated unitary 
transformation^'*. We determine the generator R of 
an appropriate canonical transformation by eliminating 
the first order term in the perturbation. To this aim, 
we impose / -I- i[R,Ho\ = 0. Introducing the nota- 
tion TZ = TZ^ + TZ^ , we decompose this constraint into 
two separate equations, which can be solved yielding: 
7^^(J-) ^ \imt^Qj-f_^r^^'>{x)dx. By this procedure^ 
we obtain the second order effective Hamiltonian for the 
magnon-conduction electron system as: 



Heff = H^and + HHe^s + - [7^^ + , + I^] 



+0 (r 



(35) 



Let us remark here that the perturbative parameter 
which actually controls this expansion, is the ratio 
Jk/Jh weighed by factors (coefficients) which depend on 
the electronic band structure and filling, as can be seen 
from the explicit expressions for and I-^ given above. 
These electronic coefficients effectively reduce the magni- 
tude of the perturbative control parameter from the raw 
value Jk/Jh, leading to a smooth convergence of the 
perturbative expansion even when | Jk / Jh \ is near or 
exceeds one. This will become clear when we present our 
results for the renormalized magnons in next section. 

Each term in the perturbation produces a correspond- 
ing term in the generator. From we obtain the "lon- 
gitudinal" generator i?^ while from we obtain the 
"transverse" generator , which are conveniently de- 



composed as: 

i?^ = ^ X! ^^rnXY -R^ = ^ RxY (36) 
X,Y=A,Bm=l,4 X.Y=AM 

The terms R^xy ^^'^ ^xy detailed in the Ap- 
pendix. 

Finally, we make a projection onto the AF fermion 
wavefunction to obtain a second-order effective Hamilto- 
nian for the magnons H'^^^. When taking the average 
(-ffe//)per„jOver the AF Fermi wavefunction we find 

{[7^^/^])p . = ([7^^,r])p . = (37) 

SO that the effective spin wave Hamihonian has the sim- 
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pier form 



^siv — (^<^ff) Fermi 



= E K (nt) + £l {nt)) + E ^^A% 



k(7 



(38) 



where 



{n^,) = [exp {£^, - /i) IkeT + l] ' (39) 
Evaluating \ (^\TZ^ ,P] + [7?.^, /^] we arrive at: 



Hl^w = E ^ + Q«) «9«9 + E («9«- 



9 + 0.qa--q 



-|-const. 



(40) 



The coefficients hQq and h}i!q are given by long and 
complicated expressions, which we detail in the Ap- 
pendix. Here we just point out that both the harmonic 
and the anharmonic parts in H'g^ have contributions 
from both longitudinal and transverse Kondo terms. 

With one last Bogolyubov transformation, we diago- 
nalize the effective Hamiltonian for the spin excitations, 
yielding: 



H 



SW 



E^^ 

9 



^^q 



qCLqUq 



(41) 



where Vlq is the renormalized antiferromagnetic spin wave 
frequency. 



III. RESULTS AND DISCUSSION 

In the previous section, we obtained a formally sim- 
ple final expression for the renormalized antiferromag- 
netic magnons in Eq. (j4ip , but which depends on a series 
of coefficients which are detailed in the Appendix. The 
long complicated expressions for these perturbatively ob- 
tained coefficients, depend on the combined effect or in- 
terplay of the different model parameters: in particular, 
some coefficients involve double and triple summmations, 
over the reduced Brillouin zone, of q-dependent filling 
factors of the conduction electron bands. Therefore, we 
have evaluated our renormalized magnon results numer- 
ically, exploring wide ranges of the different model pa- 
rameters. This has allowed us not only to assess and 
compare the influence and main effect of each of the dif- 
ferent model parameters, and show that one can reason- 
ably explain experimental magnon results^^iSi employ- 
ing parameters in the range independently shown to be 
appropriate for a phenomenological fit of specific heat 
measurements in this family of antiferromagnetic heavy 



fermionSf2£ as we will show below. Our exploration of 
wide parameter ranges has also allowed us to verify that 
the convergence radius of the perturbative series which 
determines the magnon renormalization is much wider 
than one might naively have expected. As mentioned be- 
low Eq. (pS)) . results to be shown in this section eviden- 
tiate that the actual "small parameter" controlling this 
perturbative expansion is not just the bare Jk / Jh ratio: 
this appears weighed by electronic structure and filling 
dependent coefficients which effectively reduce the con- 
trol parameter value from this ratio, leading to conver- 
gent magnon renormalization results for | Jk/Jh \> Ij 
when combined with suitable values of the other model 
parameters. 

The numerical study of our model has been done as- 
suming two simple cubic interpenetrating magnetic sub- 
lattices, for simplicity and without much loss of general- 
ity (rigorously, under this assumption our results would 
correspond to a cubic (bcc) lattice). Notice that Celna,— 
one of the compounds where inelastic neutron scattering 
(INS) on single crystals has measured the magnons we 
aim to describe, is cubic (though fee) and has a three- 
dimensional Neel-type antiferromagnetic structure. We 
evaluated magnon excitations at zero temperature, and 
assuming an underlying 3D Neel-type antiferromagnetic 
ground state of the system: our study focuses on param- 
eter sets far away from the quantum critical region of 
these systems, deep inside the antiferromagnetic phase. 
In fact, our parameters lie well inside the AF stable region 
recently determined by a DMFT -I- NRG study^^ of the 
magnetic phase diagram of the correlated Kondo-lattice 
(corresponding to the Jjj — Q case of our model: the 
addition of non- negligible AF-like RKKY coupling Jh, 
as done here, will only increase the stability of the AF 
phase). As expt. observed, in this range one might ex- 
pect Kondo-type spin fluctuations to be less relevant and 
the dispersive spin waves, object of our study, to appear 
in the AF phase. In accordance with expt. indications^ 
and for simplicity, we will further assume there is one iso- 
lated crystal field level of Ce^^ which is relevant, hosting 
a spin S = 1/2 (in fact, m = 0.5/iB is the exp. magnitude 
of the local moments in CelugSi ) ^ and concentrate on 
the numerical evaluation of the renormalized spin waves 
given by Eq. (|4ip . Hopping parameter t is taken as unit 
of energy, being W — 12t the total bare electron band- 
width. 

In the following, we will start by exhibiting and dis- 
cussing the general trends we have found in our study of 
magnon renormalization in the Heisenberg-Kondo model 
with a correlated conduction band. We shall end this 
section by focusing on the description of the measured 
magnons in antiferromagnetic heavy fermion compounds. 

Fig. Il]) depicts the first Brillouin zone (BZ) of the 
simple cubic lattice, and we include the notation for 
the special symmetry points and BZ paths on which 
the spin waves were numerically evaluated (F = O de- 
notes the zone center, A = X = 0.57r/a (1, 0, 0), Y — 
0.57r/a(0,l,0)),Z = 0.57r/a (0, 0, 1) ). Let us mention 
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FIG. 1: Simple cubic lattice: 1st Brillouin zone and special 
symmetry points. 



here that for the BZ summations we have used the spe- 
cial points BZ sampUng method by Chadi-Cohen (CC), 
at 4th order, for the simple cubic latticci^ To obtain 
dressed magnons with the correct symmetry of the lat- 
tice, we have noticed^^ it is not sufficient to use the ba- 
sic set of wave vectors of the (reduced) 1st BZ octant, 
but one needs to extend it to the (reduced) full 1st BZ 
(to properly take into account points at the frontiers be- 
tween octants): we do this applying the (48) symmetry 
operations of the Oh group to the basic set. Thus, at 
4th-order CC, taking into account all symmetry opera- 
tions we have included 5760 special symmetry points for 
each BZ summation. To achieve higher accuracy for our 
determination of the Fermi level (a delicate issue close 
to half- filling) , we have used integration over the 5th or- 
der Chadi-Cohen vectors ( that is, we used 39168 special 
symmetry points taking all symmetry operations into ac- 
count). 

In Fig. ([2]) we show the typical conduction electron 
bandstructure near half- filling, given by Eq. (|30l) . ob- 
tained with the variational approach described in sec- 
tion pi A[) . Notice that for each value of U only the 
lower band (denoted A) is filled. Chemical potential 
values obtained for the four cases shown are: /i/i = 
-0.63, 0.58, 0.49, 0.18, (slightly below the top of the low- 
est band), respectively for U/W = 0.008,0.25,0.5,0.83 [ 
i.e. U/t values of Fig. ^ ]. We obtain a direct band 
gap, determined by the AF subband energies at the cube 
diagonal BZ point: q — R (see Fig. (P) ). As expected, 
the size of the gap, as well as the energy of its centroid, 
increase with electron correlation {U) magnitude, while 
also a correlation-driven band narrowing effect is being 
described. The AF band gap value essentially depends 
on U (s), with (s) growing with both the band filling n 
and correlation U . In Table|T]these trends of the AF band 
polarization with the different parameters are recognized. 
It is also interesting to compare our AF band polariza- 
tion values (tabulated) with those reported in Fig. 2 of 
Ref. for corresponding parameters. Notwithstanding 
the different respective treatments for the band electron 
correlations, we find quite reasonable agreement where 



we could check it: in the present work we have used 
antiferromagnetic Kondo coupling values in a relatively 
narrow region around AF Jk/W ^ 10^^ — 10^''. In 
agreement with Peters and Pruschke,^'' in this parame- 
ter range we find that the local moments are almost fully 
polarized, while the corresponding band polarization ob- 
tained for U ^ and U ^ W agrees quite well with the 
corresponding values in Fig. 2 of Ref. 

12 I 1 1 1 1 1 1 1 1 

10 - - - ' - 

8 - "■■-■■-„ 




FIG. 2: Hybridized electron bands (|) along selected BZ 
paths. Parameters: S = 1/2; T = 0; i = leV\JH/t = 
0.001; Jk/Jh = l;n = 0.999. U/t = 0.1 (full lines), 3 
(dashed), 6 (dotted), 10 (dash-dotted). 

In Fig. Q we show the dependence of the AF magnon 
renormalization on the bare Jk/Jh ratio at half- filling 
( the most relevant filling for AF heavy fermion com- 
pounds): the lowest curve represents the bare magnons 
Ljq (independent of the conduction electrons). As a gen- 
eral trend, we find that the renormalization effects be- 
yond mean-field approximation reduce the spin wave fre- 
quency Vlq with respect to the MFA value Vlq, but still 
the total renormalization effect Uq Q,q is a hardening 
with respect to the bare (non-interacting) magnon fre- 
quency: ujq < ^Iq < Q.q. It is rcasouablc that due to their 
coupling to the conduction electrons, magnon excitations 
in the coupled magnetic system require a larger energy 
than when the excitation involves only the local moment 
subsystem. Thus, the dressed magnon energies increase 
with respect to the bare ones in an amount proportional 
to the increase in relative weight of the Kondo coupling 
with respect to the bare RKKY (Heisenberg) coupling, 
as Fig. ([3]) shows. Meanwhile, zJhS sets the main energy 
scale for the magnons. On the other hand, the renormal- 
ization due to only mean-field approximation effects is a 
hardening with respect to bare magnons: notice that the 
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FIG. 3: Jk/JH dependence of renormalized AF magnons: 
magnon energy along selected BZ paths. Parameters: S — 
1/2;T = 0;t = leV; Jn/t = 0.001; n = 0.999; U/t = 3. Full 
line: bare AF magnons. Jk/Jh = 0.1 (sparsely dotted line), 
1 (dashed), 5 (densely dotted), 10 (dash-dotted), 15 (double- 
dashed) . 



MFA of the longitudinal part of the Kondo interaction 
employed here, leads to a factor of zJhS+ \ Jk < s >\ 
for the "MFA-renormalized" frequency, from which the 
(longitudinal) Kondo coupling is seen to induce a magnon 
hardening: which is partially reduced when the full renor- 
malization effects are taken into account. 

Next figure ( Fig. ^ ) depicts the dependence of the 
renormalized magnons on the correlations in the conduc- 
tion band (U), at half-filling: the trend is a hardening 
of magnons when U is increased ( such hardening ap- 
pears also for the MFA magnons), and we show cases 
where U/W ranges between 0.008 (for U/t = 0.1) and 
0.83 {U/t = 10). This may be understood by the rein- 
forcement of U-dependent antiferromagnetism (spin po- 
larization) in the conduction band, which renders con- 
duction electrons less available to follow easily the spin 
excitation determined by magnon excitation in the lo- 
cal moment subsystem: an energetic cost is involved. In 
fact, the magnon gap value we obtain at g = is directly 
related to the correlation value, and to the symmetry 
breaking involved in our treatment when the itinerant an- 
tiferromagnetism of the conduction band is considered ( 
as was also reported in previous work on spin excitations 
in manganites)i^ Notice that increasing U also reduces 
q-dependent details in the renormalized magnons: this 
can be explained by analizing the indirect effect which 
U has on magnons (while Jk has also a direct effect, 
since it appears also as explicit multiplicative factor of 



0.0035 




FIG. 4: U dependence of renormalized AF magnons: magnon 
energy along selected BZ paths. Parameters: U as detailed 
here; others as in Fig. 2. Full line: bare AF magnons. U/t = 
0.1 (sparse dots: indiscernible from bare magnon curve in 
plot), 1 (dashed line), 3 (densely dotted), 6 (dash-dotted) , 
10 (double-dashed). 



the perturbatively obtained magnon corrections). U af- 
fects magnons through the modifications it induces in 
electron bandstructure, and in particular the ?7-driven 
increase of energy denominators of the perturbative co- 
efficients which determine the q-dependent renormaliza- 
tion of magnons (details in the Appendix), which explains 
the above mentioned result. The recent more refined 
DMFT-I-NRG treatment of correlations in an extended 
Kondo lattice model^i unfortunately does not allow us 
comparison, here, as their finite U results are presented 
for antiferromagnetic Jk outside the region of interest 
in our problem: their correlated AF Kondo coupling sys- 
tem is studied at much too large Kondo coupling (namely, 
Jk — 0.5W — U) for the antiferromagnetic state to re- 
main stable, being the Kondo insulator with all moments 
locally quenched the stable phase near half filling in that 
case. 

We exhibit effects of the doping on the magnon renor- 
malization in Fig. ([5]). Here the deviation of the renor- 
malized magnon energies from the bare magnon values 
is increased with the filling: at half-filling the renor- 
malization is largest, there being more conduction elec- 
trons present, which contribute to the renormalization of 
magnons by their Kondo coupling to the local moments. 
Doping away from half-filling we obtain a smooth reduc- 
tion of such renormalization effects. Thus, both filling 
and electron correlation do increase renormalization ef- 
fects, and we have already mentioned that both result in 
similar increases of spin polarization of the conduction 
band. We will shortly come back to this point, when 
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FIG. 5: Filling (n) dependence of renormalized AF magnons: 
magnon energy along selected BZ paths. Parameters: U/t — 
3; n as detailed here; others as in Fig. 2. Full line: bare 
AF magnons. n — 0.999 (double-dashed line), 0.8 (dash- 
dotted), 0.7 (dotted), 0.6 (dashed), 0.4 (sparsely dotted line: 
indiscernible from bare magnon curve in plot). 



FIG. 6: Three parameter sets leading to renormalized AF 
magnons as in Celua. (a): Jk/Jh = 100, U/t = 0.001, n= 
0.999 ( bare magnons: full line; ren. magnons: dashed); (b): 
Jk/Jh = 1, U/t= 3, n= 0.999 (dotted); (c): Jk/Jh = 1, 
U/t= 6, n= 0.7 (dash-dotted). Other parameters as Fig. 2. 



introducing our last figure. 

Let us briefly refer again to the q-dependence of the 
AF magnon renormalization we find. Some anisotropy is 
present: a larger q-dependence is noticeable along BZ di- 
agonal paths such as O — M oi O — R (see e.g. Fig. ([S])) or 
paths along the symmetry axes. While the renormaliza- 
tion effects are more pronounced at long wavelenghts: in 
particular, they are maximal at the BZ center where we 
find a spin stiffness increasing with doping, and decreas- 
ing with U or Jk/Jh- Making allowance for the quite 
different systems involved, let us mention that the renor- 
malized AF magnon behavior we obtain contrasts with 
the one recently disclosed by INS measurements in ferro- 
magnetic metallic manganites'S^ where low-q spin wave 
stiffness appears insensitive to doping, while magnons 
exhibit a doping-dependent renormalization at the BZ 
boundaries (recently suggested to be related to electronic 
correlations^) . 

At this point, let us compare our results with the few 
INS magnon measurements available for single crystals of 
antiferromagnetic heavy fermions. Comparison in more 
detail may be made only with Celua, which is cubic 
(though f.c.c.) and presents a three-dimensional AF or- 
der as we have assumed for our calculation. In Fig. ([5]) we 
show magnon results we obtained with three different sets 
of model parameters: as can be seen, all of them provid- 
ing a reasonably good description of the measured Celua 
magnons."! At the same time, these parameter sets yield 
examples of behavior of the magnon renormalization of 
the model discussed above: 



(a) the first parameter set used in Fig. 6 provides an 
example of convergence of our perturbative expansion for 
a case with Jk > Jh{— 100 Jh), at half-filling and with 
negligible U; notice a slightly larger q-dependence is ob- 
tained in this case (to be expected from the q-dependent 
expressions of the perturbative coefficients): e.g. easily 
noticeable different height of the main peaks (at symme- 
try points M and R) is visible; 

(b) with parameter set (b) of Fig. 6, again a reasonable 
qualitative description of expt. Celus magnons is ob- 
tained (considering our calculations are done for sc mag- 
netic sublattices (thus bcc lattice), while Celus in fact 
is fee). The agreement here is slightly better than with 
set (a). This shows how delicate the interplay between 
the different parameters of the model is: since, compared 
with parameter set (a), here we obtain a quite similar re- 
sult at half-filling but with Jk — Jh and U/t = 3; 

(c) finally we illustrate that we can obtain a quite sim- 
ilar description of the Celna measured magnons away 
from half filling: e.g. assuming n = 0.7 but compen- 
sating the decrease in electron energy due to doping by 
an increase in band correlation: taking U/t ~ 6 while 
keeping Jk = Jh- 

It is also worth mentioning that these sets of param- 
eters allowing description of INS results are in parame- 
ter ranges independently suggested by other authors for 
these compounds. A concrete example is the phenomeno- 
logical fit of specific heat curves for Celna , which was 
made by Lobos et al.r^ using: Jn/t = 0.0014, t = O.beV, 
n — 1 and Jk/Jh — 980; while for CeRh2Si2 they have 
used: Jn/t — 0.0034 and Jk/Jh — 430; and their data 
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extrapolation for CePd2Si2 was: J^f/t 
negligible Jk/Jh- 



0.0034 with a Physik (Univ. zu Koln) and the International Centre for 
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IV. SUMMARY 



APPENDIX. 



In the present work, we have studied spin wave ex- 
citations in heavy fermion compounds with antiferro- 
magnetic long-range order , where a strong competi- 
tion of RKKY and Kondo screening is present, as ev- 
idenced by nearly equal magnetic ordering and Kondo 
temperatures. We have described these systems using 
a microscopic model including a lattice of correlated f- 
electron orbitals (as in Ce-, U- compounds of this fam- 
ily) hybridized with a correlated conduction band, with 
the presence of competing RKKY-Heisenberg and Kondo 
magnetic couplings. 

Through a series of unitary transformations we pertur- 
batively derived a second-order effective Hamiltonian de- 
scribing the spin wave excitations, renormalized by their 
interaction with the conduction electrons. We have nu- 
merically studied the effect of the different parameters of 
this effective model on the magnon energy renormaliza- 
tion by the conduction electrons. 

We have been able to find appropriate sets of model 
parameters to describe the few existing measurements of 
magnons by inelastic neutron scattering in single crys- 
tal samples of antiferromagnetic heavy fermion Ce com- 
pounds. These parameter sets also agree with the ranges 
proposed for these compounds, through phenomenologi- 
cal fits of other experiments, like specific heat. Our re- 
sults may provide information of interest for the predic- 
tion of inelastic neutron scattering experiments in other 
compounds of this family, like CeRh2Si2, where there 
have been suggestions that the RKKY coupling should 
be stronger than Kondo effect,^''' and the only existing 
NIS measurements are of poor quality: they were made 
on polycrystals^^ many years ago, and with lower resolu- 
tion. 

As outlook towards related future work, we might 
mention the description of the experimentally reported 
magnon damping effects, and the study of the coexis- 
tence of antiferromagnetism and superconductivity in the 
context of the present model. 



In Eql32]the coefhcients M^:^^^_^ and I/^p+r-g ^i'^ a 
particular case of L^p,M^^ with A, r = ±, defined as: 

Lll = cos [4 + A^p + r (4 + ACp)] 
M^; = sin [4 + ACp + T (4 + XQ] (42) 

In EqEl the coefhcients C^y (fc, q) with X,Y ^ A,B 
are defined as: 



C+y [k, q) = SxY L++^, + (1 - 26xa) M, 



k.k-\-q 



+ (1 - Sxy) [l^U, + (1 - 2dxA) M+t+, 

Cxx 9) = ^YY 1) ^XY 9) = -'^YX 9) 

(43) 

From the longitudinal Kondo term I" one obtains the 

generator i?^ = 'Ex.y=a,b Em=i,4 R'mXY , EqEi as 
the sum of sixteen contributions. By defining the bosonic 
operator 



we have 



'.XY = ^^^^ + (1 " '^-^^) ^] (1 " ^Pl') ^ 



-ymXY Y 
^^prrj ^ra-^ p+r-qa-''-'mq 



(44) 



The operators R^xy depend on the coefficients X^-^J . 
Defining 

f^pV^ = - £l+,-q + xn + m,) 

we can write compactly 



'^pkq^ — '^xy (1 — 25xb) 



^"k,p-q+k 
pXX + - 
kpq 



Ch(^,)Ch(,9p) 



- (1 - Sxy) -pYT^Ch (z?,) Ch (z?,) (45) 

kpq 
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2XY 
pkq 



3XY 
pkq 



= (>XY (1 - 2Sxb) 



^^^k,p~q+k 
kpq 



Sh(z9,)Sh(^p) 



il-SxY)^^$y^Shi^,)Shi^,) (46) 



F, 



kpq 



= Sxy (1 - 2Sxb) 
Ltt 



^^^k,p-q+k 
pXX + + 
kpq 



Sh (i9,)Ch i^p) 



^{1- Sxy) ^$y^Sh{^,)Ch{^p) 



F, 



(47) 



kpq 
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M++ 

K^'^ = SxY (1 - 2Sxb) -p^Sh Ch {d,} 

kpq 

- (1 - Sxy) ^lh'-+ (^,) Ch (^p) (48) 



The generator i?^ resulting from the transverse 
Kondo term Eq[36l has four contributions: TZ^ — 
Y.x,Y=A,B T^XY > namely: 



T^XY = ^^\j^Y.^5xY^{l-5xY)a]y. 

kqa 

x4.n.+,.-. (Wf/at + Zf/a_ J (49) 



The coefheients and are given by: 



[Ch (^,) (fc, q) + Sh (i?,) {k, q)] 



k+g - 



(50) 



_ [Sh id,) C+y (k, q) + Ch [d,) (fc, q)] 



-"kq 



(cX_cY _ try 



(51) 



We have obtained the effective Hamiltonian in 
Eg 1381 The contribution from the perturbations 
i ([??., ^D/ermi*^^^ Written as the sum of four terms: 



1 



fermi 



q X,Y=A,B 



+ iE E {S--^+S--^){alal^ + a^a_;) 

q X,Y=A,B 

+^EM2'r+2'r)4% 

q 

+ 1 E ^ (^r + ^r) + "^a-?) (52) 



Jt^S 



= tIt E Sh (^,) (fc + g, -g) (iV^.n.,,, 



47V 
4iV 
4iV 



k(7 



J2 Wfe^^ Ch (^,) (A: + q, -q) (iV,.+,,.) 



ka 



J2 Zk^qCh id,) {k - q, q) (iVfe,-,,.) 



ka 



ka 



(53) 



q-AB 



^ J2 <f Sh (^,) C+;^ (fc + g, -g) (n^,) 



2A^ 



ka 



2N 
JlS 



J2 Wfe^f Ch (^,)C+s (fc + q, -q) {nt_,) 



ka 



f- Zt^_^C\ (^,) Cl-j, (fc - g, g) 



2N 
2N 



ka 



Zk,-qSh {i)q) C+s (fc - g, g) 



ka 



(54) 



and 
« 2iV 



5] Wfc^,^Sh i§,)C+s {k + q, -g) «,^.) 



^ E Ch i^,) (k + g, -g) «^,,) 



2iV 

' 2N 
2N 



J2 ZkA^^ i^q) c\b (k - q, q) {ntq,.) 

k 

Y Zk,-qS^ i^q) C^A ik q. q) {ntq,a) 

k 

(55) 



The S^^^ coefficients read: 



S, 



AAl 



2N 



Y <-,Ch [k ~ q, q) {Nk,-q,a) 



2N 



Y <-,Sh (^,) C+s (k - q, q) {Nk.-q,a) 
k 

(56) 



Assuming the paramagnetic band filling per site n < 1 
so that in'^„ ^ = in the ground state, one finds = 



5f ^ = 0. 



The 7^ coefficients, by defining {Nk.±q^a) 



'fc,-CT 



-^'^k±q,a) ' read: 



S, 



ABl 



E ^k:^qCh 0+-^ [k g, g) 



27V 
' 2N 



J2 <-,Sh (z9,) (fc - g, g) (n^,_ J 

k 

(57) 
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S. 



BAl 



^ <f ,Ch i^,) (fc - q, q) (n 
k 



A \ 
k — q^a/ 



2N ^ 

k 



W^^^^Sh{§,)C+^ {k-q,q) (ntq,,) 

(58) 



The coefficients S^^^ can be obtained from S^^^ by 

q q ^ 

interchanging W^Xq and Ch {-dq) respectively with Z^Xq 
and Sh i^q) . 

To write down the coefficients Vq^ and Wq^ of Eql40l 
it is convenient to introduce: 



Ch (dq 

By defining 



Sh i^q + ^p) = 6q 



M++ ^ 

_/-++ 

-M++ 

r.p— g+r 

_!_/-++ ( 
~ r.p— g+r ' 



X, 



AAl 

q,p-q+r,p 



( yABl 

(yAA3 

\ t!,P~g+r,p 



AA2 

-p,p—q+r,-q 

jAB2 



Next, by defining 



= 

"pqr 



X, 



BAZ 

q,p-q+r,p 



) £gp {''^p-q+r,a) 

_L VAB2 \ g. ( A \ 

^ '^^—p.p—q+r. — q) ^qp \ "'p—q+r,a / 

"f" '^~p,p-q+r,-q) ^ qp (j^r<y) 

&,p{nt) 



■^^r,p-q+r ^X, 



AAl 



-T++ 

r,p—q-\-r 



q,p—q-\-r,p 



x_ 



AA2 

-p,p—q+r,-q 

jAB2 



)&qp{l 



'p—q-\-r^ 



q,p-q+r,p 

r,p—q+r \ q,p—q+r,p 



yAB2 \ p, I A 

' -P.p-q+r.-q) ^IP \ "'p-q+- 



■yAAZ \ 
—P^p—g+r.-q) 



_]-+ + 

r,p-q+r 



I -yBAa 

\ q,p~q+r,p 



X 



BA3 

'P,p-q+r,- 



^qp (jT'tct) 

^gp (nt) 



and 



++ 

r.p— g+r 



AAl 

p,p-q+r. 



_T++ (y^^^ 

r.p—q+r \ —p^p—q+r. — q 



vx, 

X, 



_T++ 

r.p—q+r 



AAl 



(yA/ 

r \ q^p 



-q+r,p 
( yABi 

\ q^p-q+r-p 



AA2 

q,p-q+r,] 
BA2 

q,p—q+r,p 
AA4 

'P,p—q+r,- 

■yABi 

-P,p-q+'r,-q 



) &qp {nt) 



X 



Jqp 
Jqp 



q) ^qp 



q+' 



{''^p-q+r 



we can write 



X 



BA3 

-p,p-q+r,-q. 



(59) 



E 

pr 



a-Spq)g'^ 



(64) 



and 



-M++ ^ 

r.p— g+r 

_]-++ I 

r.p—q-\-r ' 
r.p—q-\-r 

_]-++ 

r.p~q-{'r 



(x^pj-q+r^-q + X^^^q+r,p) ^qp {^ra) in EqgOlare defined as 



Using this set of contributions, the coefficients 8^ and 



X 



-p,p-q+r. 
BAl 

-p,p-q+r,-q 



X, 



BA2 

q,p-q+r,f 



( yAAi 

\ q-p~q+'r,p 

I yABi 

\ q,p~~q+r,p 



'^~p,p-q+r,-q) '^qP " 
^AB4 



^qp (j^ra) 



I yABA \ f~ I A \ 

"T ^~p,p-q+r.-q) ^ qp \"'p-q+r,a / 

(60) 



^ ' ' X.Y=A,B 



XY 



2 ^ ' 



z+ 

q 



hV^-) (65) 



we can write 



z± 



N 



^ ] (1 3pq) C,pqr 



pr 



(61) 
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TABLE I: Conduction band AF spin polarization < s > de- 
pendence on model parameters, for the cases depicted in re- 
spective Figs. 3, 4 and 5. In first column we enter the rela- 
tive Kondo coupling magnitude (cases of Fig. 3: notice that 
Jk/W ~ Jk/ Jh X 10~^, here), while in second column is the 
corresponding band polarization < s > obtained. Similarly 
in third column we enter the electron correlation U /W cases 
from Fig. 4, and next the corresponding < s >. In 5th col- 
umn we enter filling n values of Fig. 5, and in last column 
the corresponding < s > values. 

Jk/Jh <s> U/W < s > n < s> 



0.1 0.2152 

1. 0.2153 

5. 0.2157 

10. 0.2163 

15. 0.2168 



0.008 0.0002 

0.08 0.001 

0.25 0.215 

0.5 0.406 

0.83 0.480 



0.999 0.2153 

0.8 0.1675 

0.7 0.1329 

0.6 0.0796 

0.4 0.00006 



